An approximative calculation of the fractal structure in self-similar tilings 
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Fractal structures emerge from statistical and hierarchical processes in urban development or 
network evolution. In a class of efficient and robust geographical networks, we derive the size 
distribution of layered areas, and estimate the fractal dimension by using the distribution without 
huge computations. This method can be applied to self-similar tilings based on a stochastic process. 
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I. INTRODUCTION 

Fractal nature is observed in our real-life infrastruc- 
tures of urban spatial organization and many technolog- 
ical networks. Indeed, the fractal dimensions of land- 
use and area-perimeter have been measured in real data 
of urban cities [1][2]. Similar structures with mixing of 
dense and sparse areas of nodes have been found in router 
networks [3] , air transportation networks [4] , and mobile 
communication networks [5] . These studies show that the 
spatial distribution of human activities is inhomogeneous 
and concerning with the population density. Thus, there 
commonly exists a hierarchical structure based on statis- 
tical self-similarity beyond regular mathematical models 
such as Sierpinski gasket and carpet. In other words, 
fractal behavior is not limited to objects with a regular 
morphology, but can be introduced by iterative evolu- 
tionary processes of subdivision or growing at different 
levels of scaling. 

Recently, a multi-scale quartered (MSQ) network, 
which is stochastically constructed by a self-similar tiling, 
has been proposed [6] [7] . The geographical network em- 
bedded on a planar space has several advantages [G] : the 
robustness of connectivity against failures and attacks, 
the bounded short path lengths, and the decentralized 
routing algorithm [8] in a distributed manner. Further- 
more [7] , it is more efficient (economic) with shorter link 
lengths and more suitable (tolerant) with lower load for 
avoiding traffic congestion than the state-of-the-art geo- 
metric growing networks [9-15] and the spatially prefer- 
ential attachment models with various topologies rang- 
ing from river to scale-free geographical networks. These 
properties are useful for the future self-organized design 
of wide-area wireless ad hoc networks. 

This paper investigates the size distribution of faces 
iteratively divided in the generation process of a MSQ 
network. In particular, we derive an approximative equa- 
tion for the averaging behavior of random process, and 
apply it to easily estimate the fractal dimension of such 
a self-similar tiling. 
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II. MSQ NETWORK 

For generating a MSQ network, the following process 
is repeated from an initial tiling which consists of same 
shaped faces. At each time step, a face is chosen, e.g., 
with a probability proportional to the population in the 
space of a face for the load balancing of communication 
requests, or uniformly at random. Then, as shown in Fig. 
1, four smaller faces are created in the chosen face, and a 
planar network is self-organized on a geographical space. 
Such a fractal structure is also observed in urban road 
networks [16] [17]. Note that the MSQ network includes 
a Sierpinski gasket obtained by a special selection when 
each triangle, except the central one, is hierarchically di- 
vided, however its fractal dimension log 3/ log 2 « 1.585 
differs to that in the average behavior of the following 
random selection. 




(a) triangle (b) square 




(c) triomino chair tiling (d) sphinx tiling 



FIG. 1: Basic process of the division. 



III. INFINITE STATE MARKOV CHAIN 

For simplicity, we treat the uniformly random selection 
of a face, whose case corresponds to general positions of 
nodes in the geographical network. In addition, even if 
a face in the area of high population density tend to be 
chosen for the subdivision, the assigned population to a 
node in its territory is asymptotically balanced, then the 
behavior closes to the case of random selections in course 
of time. In this setting, the subdivision process makes 
an infinite state Markov chain as illustrated in Fig. 2. 
Here, we consider a vector (ni(£), n^it), ■ ■ . , ni(t) . . .) as 
the state of Markov chain, whose element ni (t) represents 
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the number of faces at time step t on the layer I defined 
by decreasing order of size. A face on the layer I is chosen 
with the probability pi(t) = m{t)/M{t) at the next step 
t+1, then the sate is changed from (. . . , ni(t), m+i{t), . . .) 
to (. . .,m(t + 1) - l,n/+i(t + 1) + 4, . . .), where M(t) = 
Sfe n k{t) = Mo + 3t denotes the total number of faces at 
t, and Mo is the initial number. Note that pi(t) is time- 
dependent on a selected path for the decision tree from 
the top to the bottom in Fig. 2, equivalently on a se- 
lection sequence of faces. Thus, the above Markov chain 
is different from the Galton- Watson type branching pro- 
cess with a time-independent probability for generating 
offsprings [18]. 




(1,12,0,....) (2,7,4,0,....) (3,2,8,0,....) (3,3,3,4,0,....) 



(3) as shown in Fig. 3. Here, the cross, asterisk, open 
square, and closed square marks correspond to the aver- 
age of 100 states reached at each time t = 10 2 , 10 3 , 10 4 , 
and 10 5 , respectively, on the decision tree from the ini- 
tial state (3, 4, 0, . . .) in Fig. 2. The four types of dashed 
curves correspond to the distributions calculated by Eq. 
(3) until these times. The solid curve fitting with the 
closed square marks is the approximation by a normal 
distribution. Note that the averaged or calculated num- 
ber ni of faces on each layer I at t is normalized as the rel- 
ative frequency pi. This bell-shaped distribution means 
that there exist a mixing of dense and sparse areas with 
various sizes of faces whose majorities have the interme- 
diate sizes. Figure 4 shows the mean and the variance of 
the distribution pi. The upper and lower triangle marks 
denote the mean and the variance of the distribution ob- 
tained for the Markov chain. The solid and dashed lines 
denote the corresponding results for Eq. (3). In Fig. 
5, the circle marks and dotted piecewise linear line de- 
note the skewness of these distributions for the Markov 
chain and for Eq. (3), respectively. The mean and the 
variance grow as O(logt), while the skewness is almost 
constant around 0.2. Although the pi is asymptotically a 
Poisson distribution with a same value for the mean and 
the variance as mentioned in the Appendix, it is rather 
fitting with the approximation by Eq. (3) in a finite size. 



FIG. 2: Iterative division of a face chosen uniformly at ran- 
dom. 



We obtain the averaging behavior for 

Ani = m{t + l)-ni{t). 
This can be written as 

Am =4pj_i(t)-pj(t), 



(1) 



(2) 



since a face on the layer I chosen with the probability 
pi (t) is divided into four smaller ones which belong to the 
layer l + l, therefore a face on the layer / — 1 contributes 
to increase the number of faces on the layer I. For a 
large t, by noticing ni(t) = M{t)pi{t) and substituting 
M(t) =M + 3tRi3t into the right-hand side of Eq. (1), 
it is 



Am = 3(t + l) P i{t + l) -3tpi{t), 
= 3t]pi(t + l) -pi(t)\ + 3 Pl (t 



1] 



Using pi(t+ 1) » pi(t) because of t + 1 w t > 1, Eq. (2) 
is rewritten to 



Pl (t+l)-pi(t) = -j t \pi{t)-Pi-i(t)l (3) 

We numerically confirm that the average behavior of 
the Markov chain is well fitting to the solution for Eq. 
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FIG. 3: The distribution of pi. 

On the basis of the numerically obtained pi (t) , the frac- 
tal dimension df is estimated as 1.2 (a similar value is 
obtained for a chair or a sphinx tilling [19]) by using a 
cover-based method in the hierarchically coarse measure 
of doubling link length on each layer I until t = 10 8 . Fig- 
ure 6 shows an example at the state (1, 11, 4, 0, . . .). The 
number of covered areas (gray rectangles) is 

AT [2] = 2 x N[l] + 11 + 3 x (4 2 - (4 x 1 + 11)), 

where we set the number of initial coverings N[l] = 12 
for a square tiling. It depends on the primitive shape; 
N[l] = 9 for a triangle tiling, and N[l] = 8 for a chair 
or a sphinx tiling. In general for a large I, we calculate 



3 



<D 
O 
C 

!S 
to 
> 




Time step: t 

FIG. 4: Mean & Variance of the distribution of pi. 




Time step: t 

FIG. 5: Skewness of the distribution of pi. 

d } = logiV[Z]/log2 z from 

N[l] = 2 x N{1 - 1} + n, + 3 x (T, - 5,), (4) 

the total number of faces in the completely recursive di- 
vision T; = 4' , the number of holes without subdivisions 

dcf 

Si = 4xS , ;_i+n/, and Si = ni. Three terms in the right- 
hand side of Eq. (4) correspond to thick, thin solid lines, 
and dashed lines in Fig. 6. Since there is one-to-one 
correspondence between the divided four faces and the 
cross-edges in a square at each layer Z, e.g., by applying 
a clockwise mapping rule, ni represents this effect on the 
covering. However, the 2nd and 3rd terms are replaced 
by 3 x [nj/4] for a triangle tiling, and by 8 x [jij/4] for 
a chair or a sphinx tiling. We emphasis that our method 
can be performed by only using the calculated ni without 
both averaging and image processings in a box-counting 
method. 

We should remark that the continuous approximation 
with respect to time and space variables for both sides of 
Eq. (3) 
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FIG. 6: Illustration for a cover-based method. We count holes 
by the numbers for ni and the parenthesized numbers for n\ . 



gives an incorrect solution, which keeps the shape in the 
traveling wave 

Pl (t) = F{l-vlogt), 

where a function F() is determined by the initial distri- 
bution, and v = 4/3. 



IV. CONCLUSION 

Considering the stochastic network construction [6] [7] 
as a Markov chain, the size distribution of the divided 
faces has been investigated. We have shown, by numer- 
ically solving Eq. (3) as the average behavior, that the 
mean and the variance are proportional to logi, while 
the skewness is small positive. In a realistic finite size 
TV oc t < 10 8 , this approximation of pi is fitting bet- 
ter than the asymptotical Poisson solution for the in- 
teractive particle system (see the Appendix). From the 
distribution ni = Afpi, the fractal dimension is esti- 
mated as 1.2 without huge computations of averaging 
and image processings in the conventional box-counting 
method. The proposed method can be applied to other 
self-similar tilings based on a stochastic process. Fur- 
thermore, it would be interesting to discuss a self-similar 
modeling which reproduces the fractal structures of ur- 
ban road [16] [17], air-sea transportation [4], communica- 
tion [51 networks. 
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Appendix A 

Let us reformulate the subdivision process for gener- 
ating a MSQ network as a model in interacting infinite 
particle systems. Then, according to the notation [18], 
we consider the generator fin on Z + = {1,2,...}, 

tw(??)= E v(x)if(v x ' x+1 ) -mi 

where we set m = 4 and define the following functions 
on^e {0,1,2,... 



rj(x + 1) + m, y = x + 1, 

T](x) 

V{y), y^x,x + l. 



f' x+1 {y) d = { v(x) - i, y = x 



Note that the position x corresponds to the layer I of a 
chosen face in Sec. 3, therefore the number rj(x) decreases 
by 1 at x and increases by m at x due to the contribution 
at x — 1. In particular for f(rf) = 77(2), we have 



dn x {T) 
dr 



= m x h x _ 1 (t) - n x (r), x > 1, 



dn (r) 
dr 



-no(r), 



(Al) 



(A2) 



where n x (r) = E[rj T (x)] is the expectation of the num- 
ber of particles at position x and at time r, and the 
initial configuration is ?yo = (m, 0,0, . . .). We can solve 
(A1)(A2), 



n x (r) = m 1 



(x 1)! 



x > 1. 



Also, the expectation of the total number of particles is 
given by 

A/( r ) = me'" 1 - 1 ' 1 , 

Therefore, from p x (r) = ti x (t) / ! AI "(r) , we obtain the Pois- 
son distribution with a parameter mr 



Px(t) = 



(mr) 



x-l 



(x-l)\ 



x > 1. 



By the variable transformation between t and r from the 
relation N{t) = A/(0) + (m - l)t ^ N{t) = me*™" 1 ) 7 , 
the mean and the variance follow mr oc logi. This loga- 
rithmic behavior is consistent with the numerical results 
in Sec. 3. 

As another typical approach, if we consider a generat- 
ing function 



i=i 



by multiplying z l to both sides of Eq. (3) and taking the entials 
summation y\ T , we obtain 



4(1 -z) 
3t 



Here, we use ^2iPi-i(t)z l = zH^\z). The frequency 
Pi(t) is only formally calculated by the recursive diffcr- 



»® = n^zT 1 



However, it is practically unsolvable because of involvin 
with a very complicated combinatorial explosion for 
large £, equivalently for a large network size N act. 



